home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / dpptri.z / dpptri
Text File  |  1998-10-30  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPPPPPTTTTRRRRIIII((((3333FFFF))))                                                          DDDDPPPPPPPPTTTTRRRRIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPPTRI - compute the inverse of a real symmetric positive definite matrix
  10.      A using the Cholesky factorization A = U**T*U or A = L*L**T computed by
  11.      DPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPPTRI( UPLO, N, AP, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, N
  19.  
  20.          DOUBLE         PRECISION AP( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      DPPTRI computes the inverse of a real symmetric positive definite matrix
  24.      A using the Cholesky factorization A = U**T*U or A = L*L**T computed by
  25.      DPPTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  Upper triangular factor is stored in AP;
  31.              = 'L':  Lower triangular factor is stored in AP.
  32.  
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      AP      (input/output) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  37.              On entry, the triangular factor U or L from the Cholesky
  38.              factorization A = U**T*U or A = L*L**T, packed columnwise as a
  39.              linear array.  The j-th column of U or L is stored in the array
  40.              AP as follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for
  41.              1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for
  42.              j<=i<=n.
  43.  
  44.              On exit, the upper or lower triangle of the (symmetric) inverse
  45.              of A, overwriting the input factor U or L.
  46.  
  47.      INFO    (output) INTEGER
  48.              = 0:  successful exit
  49.              < 0:  if INFO = -i, the i-th argument had an illegal value
  50.              > 0:  if INFO = i, the (i,i) element of the factor U or L is
  51.              zero, and the inverse could not be computed.
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.